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We consider a comprehensive set of simplified models that contribute to final 
states with top and bottom quarks at the LHC. These simplified models are used to 
^-H create minimal search strategies that ensure optimal coverage of new heavy flavor 

physics involving the pair production of color octets and triplets. We provide a 
^ ^ set of benchmarks that are representative of model space, which can be used by 

O experimentalists to perform their own optimization of search strategies. For data sets 

larger than 1 fb~^, same-sign dilepton and 36 search regions become very powerful. 
OO Expected sensitivities from existing and optimized searches are given. 
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I. INTRODUCTION 

Data from the LHC has begun to inform our understanding of electroweak symmetry 
breaking (EWSB). Since naturalness of EWSB typically requires the existence of a top 
partner, many promising new physics scenarios produce final states rich in top and/or bottom 
quarks (heavy fiavors), see e.g. lU [2]. In many cases this is accompanied by a stable neutral 
particle in the final state that could be the dark matter of the Universe. The 7 TeV LHC 
run has the potential to discover new colored states with masses well into the TeV scale 
that give a new contribution to events with heavy fiavor content. In this paper, we study a 
comprehensive set of simplified models to aid the 7 TeV LHC search for new colored particles 
that decay to top (t) and/or bottom (6) quarks and to a new stable particle that appears as 
missing energy. 

A generic difiiculty in LHC new physics searches is that the masses of the new states 
produced are not known a priori. This means that the kinematics of the final state signatures 
can vary drastically, depending on the mass spectrum of the particles produced. For heavy 
fiavor final states, an additional challenge in developing a comprehensive search strategy is 
that top quarks lead to different signatures than bottom quarks. For example, a 4t final 
state is quite distinct from a Ab final state: the former produces a large multiplicity of visible 
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particles, thus increasing (decreasing) the energy (missing energy) of the event, while the 
latter has a smaller multiplicity of visible states and a large missing energy from the stable 
particle at the bottom of the decay chain. Given finite resources, it is important to design 
a limited set of search regions that nevertheless retain close to optimal discovery power for 
all possibihties. 

An effective but minimal set of search regions can be developed using "simplified models" 
[3H6], which are effective field theories with a minimal particle content. Only the physically 
relevant parameters of a full model are kept, namely the masses of the new states, the 
production cross section, and the branching fractions into the available modes. Different 
simplified model parameter values can lead to distinct signatures and thus require disparate 
search strategies for discovery. Simplified models allow one to efficiently explore all pos- 
sibilities and cast a wide net for physics beyond the Standard Model (SM), a necessity in 
searches for new physics due to the LHC inverse problem [7] . 

This article will consider 12 simplified models in which either color octets, g, or triplets, 
q, are pair produced and decay to a combination of t- and/or 6-quarks and a stable particle 
X^, with possibly some light fiavor jets. Early LHC searches are already looking for these 
types of spectra, which offer an optimistic path for discovery because of the highly efficient 
6-quark tagging at the ATLAS and CMS experiments. 

From the theoretical side, as alluded to above, these signatures are well-motivated because 
the mechanism stabilizing the mass of the Higgs boson typically relies on a new partner 
for the top. Moreover, naturalness arguments favor light third generation partners while 
constraints from fiavor changing neutral currents and CP violation typically favor much 
heavier light fiavor partners. A representative example is supersymmetry, in which the 
gluino is the color octet that can be pair-produced and will decay predominantly to t and 
b quarks through the light stops and sbottoms if these are much lighter than the first and 
second generation squarks. One of the advantages of simplified models is that the results can 
be recast into different theories. For example, universal extra dimensions (UED) theories 
have a similar particle content to supersymmetry, but different spin quantum numbers. This 
last difference can be accounted for by a re-scaling of the production cross section by the 
appropriate degrees of freedom, as long as spin-correlation effects are subdominant [S]. 

Within each of the 12 simplified models, we consider a discrete but broad set of masses 
for the g (or q) and the giving us a total of 2762 models. We first find the 'best' 
search strategy for each of the 2762 models, which consists of the combination of cuts 
that optimizes the square-root signal over background. We consider 8064 combinations of 
cuts on the number of jets, b-jets, leptons, missing transverse energy (-^t), and total event 
transverse energy {Hj)- We then find a set of 'good' search strategies, namely a set of 
counting experiments that will discover all 2762 models with great efficiency. The 'efficacy' 
£^ of a particular search region is the ratio of the quantity of data needed to discover the 
model with this search region to the amount of data needed using the 'best' cut. Using a 
genetic algorithm, we identify a small set of search regions ((9(10 — 20)) that cover all 2762 
models with £ < 2. Noteworthy is that 36 and same-sign dilepton search regions (SSDL) 
are particularly sensitive to these models, especially with larger data sets. 36 search regions 
can become an effective way to discover a broad class of heavy flavor topologies from g 
pair production, while SSDL signal regions can help in multi-top topologies where the mass 
spectrum of new states is compressed. 

As theorists, we do not have the best tools for simulating the LHC detectors, so the 
set of search regions found in this study may not be optimal. Therefore the real goal of 
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this article is to suggest some general lessons and to provide a methodology for optimizing 
searches over a large parameter space of new physics models. Moreover, we identify a set 
of benchmarks consisting of (9(100) simplified models that are chosen from the full set of 
2762. These benchmarks can be used by LHC experimentalists to optimize their search 
strategies and still retain excellent sensitivity to the full space of models. In other words, 
these benchmarks taken together are representative of all 2762 models, and optimizing search 
strategies on the benchmarks is to a good approximation equivalent to optimizing on all 2762 
models. This approach significantly cuts down the work that needs to be done by the LHC 
experimentalists in picking sensible benchmark points, and also minimizes the possibility 
that there are regions in parameter space that might be missed. An important caveat here 
is that the optimization is done over a particular set of cuts (see above), and we cannot 
assure full model space coverage if different cuts are used. However, the benchmarks chosen 
offer vastly different final states. The fact that (9(10) search regions are needed to cover 
the space of 60 benchmarks is a statement that each benchmark can give quite different 
kinematics from the rest. 

The organization of this paper is as follows. Sec. [TT] introduces the simplified spectra used 
in this study. The signal and background calculation is discussed in Sec. HI The reach of 
the LHC's latest analyses on the simplified models is estimated in Sec. IV A minimal set of 
search regions is found and presented in Sec. |Vj Sec. VI summarizes our results. Appendix 
[A] and [B] show the expected cross section sensitivities for our simplified models from current 
ATLAS searches and optimized searches, respectively. Appendix [C]presents the benchmarks. 



II. SIMPLIFIED MODELS 



A variety of search regions are necessary to comprehensively cover new physics scenarios 
that produce heavy flavor final states. One reason for this is that different numbers of t 
versus b quark final states can significantly alter the optimal search strategy. We introduce 
12 simplified models, which we divide into two broad classes: (i) gluino-like g, where the 
g are color octets that are pair-produced, and each decays to two third generation quarks 
and a light electroweakino-like state ('LSP'); and (ii) squark-like q, where the q are color 
triplets that are pair-produced, and each decays to one third generation quark and a light 
electroweakino-like state. In the g simplified models, we always assume the g decay through 
heavy off-shell squark-like particles, which themselves are inaccessible at the LHC. Similarly, 
in the q simplified models, the g are assumed to be much heavier than the q. We also include 
light flavor decays of the g for completeness. 



A. Gluino-like Models 

We divide the g simplified models into two further categories, characterized by the nature 
of the LSP. In the first category (referred to as topologies), the LSP consists of a single 
neutral state, x^- There are abundant examples of such theories, such as supersymmetric 
theories where the lightest supersymmetric particle is a bino or singlino. 

In the second category (referred to as Q^^ topologies), there is a charged state nearly 
degenerate with the neutral state The transition ~^ X^y which occurs through an 
off-shell W^, is effectively invisible at the LHC. In all Monte Carlo (MC) calculations, we 
set m^± —m^o = 10 GeV, which is an unobservable mass splitting in jets and missing energy 
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FIG. 1: Gluino-like simplified models used in this study. Qf^ models have only a light 'bino'- 
like state, while Qy^ models have an additional charged state. The charged state is always nearly 
degenerate with a neutral state, so that there are not any additional visible leptons from its decay. 



analyses. In supersymmetric theories, this spectrum frequently occurs if the LSP is a Wino 
and the Higgsinos are heavy or if the LSP is a Higgsino and the Wino and Bino are heavy. 
More generally, a degenerate ^ind are well-motivated whenever the LSP is part of an 
SU{2)l multiplet. 

We now map each possible decay of the gg pair to a simplified model. 

1. Topologies 

The g in the Qj^ spectra can decay in various ways, as illustrated in the left panel of 
Fig. [T] The modes are 

T: (la) 

B: -g^x'hh, (lb) 

J: g^X^JJ- (Ic) 
Since the g are always pair produced, the following six decay topologies are possible: 

r'TT r'BB r'TB r'TJ r'BJ nJJ fn\ 
^B' ^S' ^S' ^B' ^B' ^B- (2j 

A thorough study of the Q'j^ topology was performed in [B]; our study of light flavor jets 
here will be cursory and only done to allow comparisons with the heavy flavor searches. 

The and Q'j^ will be referred to as pure topologies, and in principle these decay 

topologies could be the only signal of new physics. This is in contrast to the hybrid topologies 
Qj^, , which each necessarily will be accompanied by another topology. For example, 
if the topology is present, the Q"^ and Q^J topologies must also be present in the data. 
Although hybrid topologies never appear as the only topology, they may have the largest 
branching ratio and can frequently be the most visible. For instance, if Br(^ — > T) = 10% 
and Br(^ — )■ J) = 90%, then the topologies appear in the ratio 

• • = 1% • • 81%- (3) 
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However, the appearance of final state top quarks can make mucli more visible than Q^^^, 
compensating for the smaller rate. We will develop optimized strategies for both pure and 
hybrid topologies, regardless of the branching ratios into particular final states. 

2. Qyi^ Topologies 
Models that have a Q^^ spectrum have a new decay mode 

M : g X' x^ibW^* (4a) 

M : g^x^bt^x^biW-\ (4b) 

where in all cases the boson is so far off-shell that it is effectively invisible (as mentioned 
above, for definiteness, we choose it to have 10 GeV of energy). Note that > 100 GeV 
from chargino searches at LEP-II [S] , although in plots we will allow for much smaller . 
We always take M and M in equal admixtures and drop the distinction between the two decays, 
since they look very similar in detectors. There are still the direct decays of g to which 
can give rise to T, B and J decay modes. Constructing the ten possible decay topologies 
from these four decay modes, six of them are identical to decay topologies. The four new 
decay topologies are 

y^, t/^, yy^, y^- (Sj 

Eqs. ([2]) and ([s]) thus give the ten possible decay topologies in g simplified models, when 
cascade decays are either suppressed or inaccessible, that we will study. In supersymmetry, 
we could see various specific admixtures of these topologies, depending on the spectrum of 
squarks and the identity of the LSP. 

The decay width of a gluino into a fermion species is T~^^^- oc /m'^, where Yg is the 
quark hypercharge. Of course, due to the large t-quark mass, phase space considerations are 
important and can significantly modify the branching ratios. Let us consider the branch- 
ing ratios in several examples of plausible supersymmetric topologies, ignoring phase space 
factors: 

• Bino-like LSP, lighter right handed third generation squarks. Assuming that m^c = 
m^c, this gives a ratio of g decays of T : B = 4 : 1 resulting in 

gl^ : : gf = 64% : 32% : 4%. (6) 

• Wino-like LSP with TT decay mode kinematically accessible 

g^ : gf : gf ■. g"^ ■. g^^ ■. g^^ = 2.8% : 5.6% : 2.8% : 44% : 22% : 22% (7) 

• Wino-like LSP with TT decay mode kinematically inaccessible 

gf ■.g^^:g'^ = 11% : 44% : 44% (8) 
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FIG. 2: Squark-like simplified models used in this study. 



B. Heavy Flavor Squark-like Models 



The second class of simphfied models consists of a single color-triplet squark-like state, 
either t or 6, which are pair-produced. There are again two categories characterized by the 
nature of the LSP. A single neutral stable particle, allows the decays 

t: 

b: (9) 

giving rise to simplified models that we call 7^ and S^, respectively, with topologies T| and 

Supersymmetric theories with comparable t and h masses can give 0{1) admixtures of 
these two processes, but we will not consider these scenarios here. 

In 7^ and Bt^ simplified models, a chargino-like state nearly degenerate with the x° 
LSP allows the decays 

t: h^x't^x'^tW-^ 

b: t^x'^b^x^bW^*, (10) 

where we assume the decay products are soft and unimportant, and we therefore again 
label the decays as t and b as in Eq. Various decay topologies are possible, but they 
are again effectively the same as 7^ and B\ considered above. 



III. BACKGROUNDS AND SIGNAL SIMULATION 

The dominant SM backgrounds to jets and I^t signatures are ti + jets, + jets, 
+ jets, t + jets, VV+ jets, tV jets and QCD, where V = Z^,W^. The matrix 
elements for parton level events were computed in MadGraph 4.4.32 [10] with CTEQ6L1 
parton distribution functions [TT]. Variable renormalization and factorization scales are set 
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FIG. 3: Feynman Diagrams for t + nj and t + W + nj. Notice the need to subtract the on-shell 
top contribution from tWb sample. 



to the transverse energy of the event [12]. The SM parton level processes generated are 

pp^W^ + nj l<n<3 (11) 

pp ^ + nj 1 < n < 3 

pp ^ ti + nj < n < 2 

pp — )■ VV + nj < n <2 

pp — > t/i + nj 1 < n <3 

pp ^t/t + W^ + nj 0<n<2 



where j G {g, u, u, d, d, c, c, s, s, b, b} . 

The SM contribution to distributions peaks at low energies, whereas many new physics 
signatures produce events with large ^t- Therefore it is important to have sufficient MC 
statistics on the tail of the distribution. To achieve sufficient statistics, different samples 
are generated for each SM process, where each sample has the pt of the massive particle 
lying in a given interval. For instance, Z° + jets parton level events were divided into three 
samples with 

GeV < Pt,zo < 200 GeV 
200 GeV < pt,zo < 300 GeV 
300 GeV < pt,zo ■ 

In the case of two massive particles produced, such as tt, the samples are divided by the 
larger px of either massive particle in the event. 

Contributions from QCD to jets and can come from either detector effects and jet 
energy mis- measurement, or neutrinos appearing in the decay of heavy flavor hadrons. To 
estimate the QCD contribution to jets and signatures, the following subprocesses were 
generated 

pp — )■ nj' 2 < n < 4 
pp->bb + nj <n<2. (12) 

Here j' refer to light flavor jets and gluons only. To achieve sufficient statistics, the QCD 
and bb backgrounds were subdivided in exclusive bins delimited by the pt of the leading jet. 
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The signals, g and q pair-production, were generated in association with up to two jets 
at parton level. 



PP-^ 99 + nj 
pp ^ qq + nj 



<n<2 
<n<2 



(13) 



The effects of including additional radiation in signal processes have been documented in 
several studies [H El [131 HI]- We generated MC for 2762 different points in the rrig 
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planes. 



For both signal and backgrounds, the showering, hadronization, particle decays, and 
matching of parton showers to matrix element partons are done in PYTHIA 6.4 [T3]. We 
use the MLM parton shower /matrix element matching scheme with a shower- scheme 
introduced in [13]. A fixed 5-flavor matching scheme is used. The matrix elements better 
describe hard radiation, while the parton shower generates softer radiation that fills out jets 
|14j . The matching scales used here are: 



Sample 


Q Match 


tt + jets 


100 GeV 


\/ + jets 


40 GeV 


t + jets 


100 GeV 


tV + iets 


100 GeV 


QCD 


50 GeV 


99 + jets 


100 GeV 


qq + jets 


100 GeV 



(14) 



Hard jets beyond the multiplicities listed in Eq. (11) must be generated by the parton 
shower. In particular, for + jets and Z° -|- jets, the fourth jet and beyond are generated 
through the parton shower. This approximation has been validated by several studies [H]. 
As a cross check, for the Z° -|- jets background generated for our study, the discrepancy in 
the inclusive rate for four jets and from matching up to three jets versus matching up to 
four jets is C(15%), assuming a selection requirement of 50 GeV on the fourth leading jet. 
For consistency, in the analysis only samples with up to 2 — )■ 4 partons are used. 

The matching scheme is particularly important for the signal when the spectrum of the 
new states is compressed. In this signature, the final state jets from the signal can be soft 
and mimic QCD events. The signal can be enhanced by requiring hard jets, which will come 
from ISR or FSR. Fig. |4] considers the pt spectrum of the leading jet and the jet multiplicity 
from a 400 GeV g that decays to a degenerate 390 GeV x° in the topology. Significant 
differences are seen, generally with the unmatched sample underestimating the tail of the 
distributions. 

Next-to-leading-order (NLO) corrections alter the predictions of both signal and back- 
ground. With parton shower /matrix element matching, the shapes of differential distribu- 
tions are accurately described by tree level predictions. The largest corrections are to the 
inclusive production cross section and can be absorbed into fC-factors. The leading order 
cross sections of the signal are normalized to the NLO cross sections calculated in Prospino 
2.0 [16]. Fig. [5] shows the NLO cross sections for both g and q pair production. The 
ti + X,W'^ + jets, and + jets leading order production cross sections are scaled to the 
NLO ones from Hj 
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FIG. 4: The effect of matching is shown for a signal with a 5 at 400 GeV and a x at 390 GeV in 
the Q^^ simphfied model. The left plot shows the pT spectrum of the leading jet. The right plot 
shows the effect on the jet multiplicity. 
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FIG. 5: NLO cross sections for gluinos and third generation squarks, with other squarks decoupled. 



PGS 4 is used as a detector simulator [18j. We use the PGS 4 ATLAS card, which has 
been shown to reproduce results to C(20%) accuracy. 

One of the drawbacks of PGS 4 is that it uses a cone jet algorithm with AR = 0.7. This 
is an infrared unsafe jet algorithm, but better represents the anti-/cT algorithms used by the 
experiments than the kx algorithm. The SM backgrounds change by at most (9(10%) when 
varying the cone size to AR = 0.4. 

The signal offers a more varied testing ground for the effects of changing the jet algorithm. 
Two competing effects are found. The first is that there is more out-of-cone energy for 
smaller cones, resulting in less energetic jets. The second effect is that smaller cone jet 
algorithms find more jets. The dependence of the kinematic cut efficiencies on AR varies 
with mass splitting between the g and x- compressed spectra, when the pr of the jets 
is reduced, the efficiencies for the smaller cone size decrease because jets fall below the 
minimum jet px requirement. For widely spaced spectra, where jets are energetic, more 
jets are found with a smaller cone size and the efficiency to have multiple jets passing the 
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minimum jet requirement increases. Altogether, the efficiencies differ by at most (9(20%) 
and are consistent with other studies [19j. This effect is not included, and we simply use a 
fixed AR = 0.7. 

We have modified the publicly available PGS 4 code to include a different algorithm for 
identifying b- and charm-jets. The algorithm now matches 5-mesons at the MC truth level 
in the PYTHIA output to the jet found by PGS 4 that lies closest to it in AR. In this way, 
all 5-mesons will be matched to some jet. This jet is then preliminarily identified as a 6-jet 
if it lies within AR < 0.4 of the 5-meson. Note that different i?-mesons (assuming there 
is more than one in the PYTHIA event record) can be matched to the same jet, so that the 
number of 6-jets in an event can be lower than the number of 6-quarks produced in the 
event at the generator level. We checked the efficiency for finding 6-jets in this way in two 
different samples. For a sample of pair-produced sbottoms with a mass of 400 GeV, where 
each sbottom decays to a bottom quark and a neutralino, we find an efficiency of ~ 86% 
(89%, 90%) when requiring both 6-jets to have px > 20 GeV (10 GeV, GeV). Note that 
the efficiency increases by 1% if we increase AR between the jet and the 5- meson to 0.7 
from 0.4. The efficiency of the original PGS 4 algorithm is 83-84%. For a sample of W^bb 
events, we find an efficiency of 75% (78%, 79%) when requiring both 6-jets to have pt > 20 
GeV (10 GeV, GeV). 

After the 6-jets are identified, we identify c-jets in a similar way, except now D-mesons 
are matched to jets that are not already 6-jets. 

Having identified the 6- and c-jets at the MC truth level, we randomly turn a 6-jet into 
an ordinary jet without a 6-tag with a probability given by the px and ?7-dependent 6-tag 
efficiencies found in [20|. We also turn ordinary jets and c-jets into 6-jets with a pT-dependent 
mistag rate given in [20] (for the mistag rate, we average over the small r] dependence). 

With these modifications to the PGS 4 code, the main standard model backgrounds were 
calculated. Fig. |6] shows the pt spectrum of the leading, second, and third 6-jet in the event 
for each of the SM backgrounds that were calculated. Fig. [7] shows the distributions for 
backgrounds after requiring > 1 and Ub > 2, respectively. Furthermore, Fig. |8] shows 
a comparison between the main SM backgrounds generated for this study in PGS vs the 
simulated backgrounds used by ATLAS in [2I| • Reasonable agreement between distributions 
is found for rif, > 1. For > 2 the distributions obtained from PGS over-predict the number 
of events in the signal region. ATLAS currently does not use a separate jet energy calibration 
for 6-jets which are systematically lower reconstructed at lower pt than their true value. To 
model this effect, the pt requirement on the 6-jets is increased from the quoted value by 
ATLAS of 50 GeV to 60 GeV and gives a better agreement between our backgrounds and 
those from ATLAS. 

IV. EXPECTED LIMITS FROM EXISTING LHC SEARCHES 

The LHC experiments have performed analyses in the jets and channel using nearly 1 
fb~^ of data from the 2010 and 2011 runs. It is important to study and compare the search 
regions in these studies to our optimized search strategies in Sec. |Vj We consider several 
ATLAS studies, three with 35 pb~^ of 2010 data that remain relevant for constraining the 
low-mass regions, and three that include more data: 

1. A study [22] with 35 pb^^ of 2 or 3^ light flavor jets and > 100 GeV, with various 
meff and jet px cuts, including cuts on the ratio of ^x/^eS 
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FIG. 6: The px of the 1st, 2nd, and 3rd hardest b-jet in MC background events with at least that 
many &-jets, with a total luminosity of 1 fb~^, after requiring 100 GeV. 
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FIG. 7: of MC background events with greater than or equal to 1 or 2 6-jets, with a total 
luminosity of 1 fb~^. 



2. A study [23] with 35 pb ^ demanding at least one fe-jet with or 1+ lepton, and 

> 80 or 100 GeV, respectively 

3. A study [21] with 35 pb~^ using opposite-sign or same-sign dileptons, with a missing 
energy cut ^t> 150 or 100 GeV, respectively 

4. An updated study [25] of jets, leptons, and missing energy with 165 pb~^ requiring 
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FIG. 8: Comparison of the backgrounds generated for this study against the ATLAS backgrounds 
from [21j. The top (bottom) panel shows the f^T (-^eff) distributions in the n^y > 1 and Uf, > 2 
channels, respectively. The backgrounds generated in this study with PGS are shown in solid black, 
the ATLAS backgrounds in dashed black, and the data in solid red. The figures are shown assuming 
C = 0.83 fb~^ of integrated luminosity. 



> 0.25mcff and > 500 GeV 

5. An ATLAS study |21] with 0.83 fb^^ requiring ^r/m^s > 0.25 and > 130 GeV, 
defining four signal regions with 1 or 2 6-jets and rrics > 500 or 700 GeV 

6. An updated study ^26j with 1.03 fb^^ demanding rrics > 1000 GeV and ^T/f^eS > 0.25 
with four signal regions with various jet pr cuts 

There are also many similar relevant CMS analyses looking for jets and missing energy with 
or without leptons and photons, analyses using MT2 and ax, and also a CMS analysis 
using same-sign dileptons (e.g. [271 - EU] ). While this paper was in preparation several new 
interesting LHC searches have emerged [30]. We will not consider these here. 



In Appendix |Aj, in Fig. [14| we show the expected 95 C.L. limits for £ = 1 fb~ from 
the above ATLAS studies on all simplified models from Sec. [TT] for different choices of the 
production cross section. Note that we use our own MC for the background modeling. In 
the next section, we find a set of minimal search regions that are needed to cover the space of 
simplified models, and estimate the 7 TeV LHC's sensitivity to the simplified models using 
this minimal set of search regions. 
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V. OPTIMAL SEARCH REGIONS 

The simplified models introduced in Sec. [IT] can be used to develop broad search strategies 
that cover the model space. Despite the reduction in the number of relevant free parameters 
in simplified models compared to complete theories, a multi-signal-region strategy is needed 
to efficiently cover all kinematic possibilities. In this section, we find a minimal number of 
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signal regions necessary to cover the entire space of simplified models. Then in Sec. VB 
we will propose a set of benchmark models that span the full parameter space of simplified 
models, in the sense that a search strategy that is sensitive to all benchmarks will also be 
sensitive to all of the simplified models. 

The terminology used throughout the rest of the article is introduced in what follows. 
We assign all events to a signal region defined by the number of jets (iVjet), 6-jets (iV^jet), 
and leptons (iV^), as well as the missing transverse (^t) and visible energies {Ht): 

{N^,„Ne,N^iet,:^T,HT). (15) 

In addition to a cut of 120 GeV on the transverse momentum [p^) of the hardest jet, we 
consider the following set of cuts: 

{2+,3+,4+} 
{0+,l+,2+,3+} 

{0, 1+, 2+, 3+, SSDL+, OSDL+} 
{0,50, 100,..., 500} GeV 
{200, 300,..., 1200} GeV. 

Here labels without a "+" are exclusive cuts, e.g. exactly zero isolated leptons are required 
for Ni = 0; the superscript indicates that the search regions are inclusive, e.g. N^et = 2"*" 
requires two or more jets; SSDL refers to "same-sign di-lepton" and OSDL refers to "opposite 
sign di-lepton"; leptons are electrons and muons (not taus, which are treated as jets) and 
are required to be isolated; all jets (including 6-jets) beyond the hardest one are required to 
have pt > 50 GeV and |?7j| < 2.5. This gives rise to a set of 8,712 possible search regions. 
Not all search regions are physical due to the overlap between the Ht cut and other pt 
and/or requirements and this reduces the number of search regions to 8064. Note that 
requiring even higher jet multiplicities may be useful [3Tj, but we do not consider this here. 

A given signal region or cut, Cj, will yield an expected limit on the cross section times 
branching ratio, a x B, for a given simplified model at the 95 % C.L. given by 

(a X B), = ,^^^,]' . (16) 

where e(M)j is the efficiency of Cj on the model point M. A{B)i is the allowed number of 
events in the signal at the 95 % C.L. if B background events are expected and in fact fit the 
data. We take 



A{B) = 2 X ./Stat(5)2 + (e(nbjct)syst5)2, (17) 



where Stat(i?) is the Poisson limit on B. We also include a systematic error, e(nbjct)syst5 as 
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a function of N^et in the signal region ^: 

e(06)syst = 20% e(16)syst = 20% e(26),yst = 40% e(36)syst = 60%. (18) 

The optimal cross section limit on a model M is given by 

(a X B)opt = {min((a x B),) : t e {1, N,^ts}} , (19) 

where the number of search regions is Ai'cuts = 8064. It is natural to quantify the "goodness" 
of a cut Ci by the amount of LHC data needed to make a discovery or exclusion using that 
cut. For this purpose, we introduce the efficacy of a cut 

m) = (20) 

[a X B)opt 

In words, this is just the ratio of the necessary production cross section for discovery using 
a cut Ci divided by the cross section necessary for discovery using the optimal set of cuts. 
An efficacy of 1 is 'perfect', and otherwise smaller efficacies are better. Thus the best search 
strategy for all model points M will be a combination of cuts {Ci} such that £ is close to 
one for every model using at least one of the cuts in the search strategy. 

It is interesting to compute the expected efficacies of the search regions used by public 



LHC analyses (see Sec. IV) for the heavy flavor simplifled models. Fig. |9] shows the efficacies 
for the Q"^^ (top panel) and (bottom panel) topologies for £ = 1 fb~^, 5 fb~^, and 15 
fb~^ (flrst, second, and third column, respectively). For the current low-luminosity searches, 
the LHC analyses have very good efficacies, and there are even a few small isolated regions 
where the LHC searches have been slightly better than the search regions considered here, 
due to the looser triggers allowed by lower luminosity analyses. At higher luminosities, it is 
of course possible to greatly improve the search strategies. 



A. Optimizing Search Strategies 

It is not feasible, nor necessary, to look at all 8064 search regions to have a reasonable 
efficacy over the entire set of simplifled models. One of our principle goals is to flnd a minimal 
comprehensive set of signal regions that cover the space of simplifled models spanned by the 
topologies from Sec. [ll] Each solution to this problem will be a set of cuts {Ci}, where 
at least one Ci in the set gives S{Ci) < Sent, for some "critical" efficacy Sent, for every 
point M in the space of simplifled models. Different solutions will be found for benchmark 
luminosities C = 1 fb~^, 5 fb~^, or 15 fb~^. 

A genetic algorithm was used to construct the minimal set of search regions to cover 
the entire space of simplifled models. The conflguration space for the genetic algorithm is a 
binary string of 8712 bits that signify whether a particular search region is used or not (recall 
that there are 8712 search regions before imposing constraints); thus the conflguration space 
is 2^^^^ states. A set of random cuts were created by turning on a weighted selection of 40 
of the 8712 search regions. The weight of each search region is proportional to the number 



The choices of systematic errors used were made fohowing private communication with the SLAC ATLAS 
group; however, these are also consistent with |21j . 
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FIG. 9: Expected efficacies of the ATLAS search regions on the G"^^ (top panel) and G^ (bottom 
panel) simplified models, for integrated luminosities C = 1 fb^^, 5 fb^^, and 15 fb~^ (first, second, 
and third column, respectively). Our goal is to find a set of search regions that have efficacies 
close to 1 for all simplified models and for each luminosity. While current ATLAS searches give 
excellent coverage, other strategies will be required for higher luminosities data sets, especially for 
events with many final state 6-jets and tops. 



of models to which the search region is sensitive (i.e. the number of models for which the 
cut Ci gives £{Ci) < f^crit)- The sets of search regions are evaluated to see how many models 
they cover within the desired efficacy, and a "fitness" is assigned to them with the formula 

nC,M)^j,^-^y^^, (21) 

where M is the number of models covered, C is the number of search regions, and M^ax is 
the total number of models. This fitness function strongly penalizes search strategies that 
do not cover all models, followed by a penalty for having too many search regions. 

After evaluating the fitness of the search strategies, the least fit 50% are removed. Pairs 
of fit search strategies are then selected and a new search strategy is created by taking a 
randomly determined fraction of each search strategy's used search regions. For instance, 
if the two selected search strategies had A^^i and N2 search regions, then a uniform random 
number on the unit line segment, x, would determine that xNi search regions would be 
taken from the first search strategy and {1 — x)N2 would be taken from the second search 
strategy. So if Ni = 20 and A^2 = 30 and x = 0.20, 4 search regions would be taken 
from the first search strategy and 24 would be taken from the second. If duplicate signal 
regions are selected, the duplicate is removed, reducing the number of search regions. After 
creating a new search strategy, the search is mutated to guarantee that the population of 
search strategies had sufficient diversity. Each used search region has a probability of being 
changed to another random search region. We use 6% for this probability known as the 
"mutation rate". Thus with the 28 search regions in the example, 1.5 changes would be 
made on average. 

If after ten consecutive generations no progress has been made, i.e. if no solution has been 
found that covers the entire model space, then a solution is manually created by forcing every 
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model to be covered by some search region. This can be done by increasing the number of 
search regions in the search strategies until full coverage is achieved. Finally, if every model 
is covered and no further progress is achieved for seven generations, search strategies are 
scoured to see if any search region can be removed without reducing coverage. Either way, 
the genetic algorithm is restarted. If no progress in reducing the number of search regions 
in a search strategy has been made in twenty generations, the program ends. 

Typically, the program terminated after 50 to 70 generations, and 20 to 300 distinct 
optimized search strategies were found each time. While the termination of the program does 
not guarantee that the optimal solution has been found, re-running the program multiple 
times usually results in the same number of required search regions. The distinct search 
strategies typically have similar features even if they differ slightly in detail. 

The program can easily identify optimized search strategies that cover all models with 
an efficacy less than the critical efficacy, and we found that Ecrit = 1-75 results in just 
a small number (< 10) of signal regions that can cover all models. We have chosen to 
present our benchmark optimized search regions, which will be discussed in more detail in 
the next section. Table [T] presents the search regions found after optimizing over all 8712 
search regions that together cover all benchmark models with an efficacy of 1.75 or better 
(i.e. £crit = 1.75), for £ = 1 fb-\ 5 fb-^ or 15 fb-^. 

The optimized search strategy from Table |T] has three distinct sets of search regions. 
First, it achieves sensitivity to generic light flavor jets and missing energy signals with a 
small number of regions that involve significant and large Ht cuts. Second, to uncover 
heavy flavor physics, there are several search regions involving 6-jets with more modest cuts 
on Ht- This second category of search regions is not that different from standard signal 
regions used by the LHC experiments thus far in searches for heavy flavor models. Third, 
a set of non-standard cuts are found. Since many of the simplifled models produce 46 jets 
in every event, we flnd that 36-jet search regions with a very modest cut and a minimal 
Ht cut can be very effective (see also [2]). Furthermore, the signal region involving same- 
sign dileptons in events with 2 or more 6-jets achieves great sensitivity to multi-top events, 
particularly for points in model space where ~ 2mt + m^o. 



B. Benchmark Models and Re-optimized Search Strategies 

The optimized set of signal regions discussed above cannot be obtained without sam- 
pling over a very large model space. There may be practical limitations in doing such flne 
sampling, especially for the experimental groups which must use full detector simulation in 
their analyses. In this section, we flnd a set of models, or benchmarks, so that if a search 
strategy that covers this set of models with a given efficacy, ^crit, is found, it will also cover 
the entire space of simplified models. An important caveat here is that this assumes that 



they optimize over the set of cuts in Eq. (15). Additional cuts, or a different set of cuts, 
may require slightly different benchmarks. 

To make the set of benchmarks as intuitive as possible, we began with five benchmarks per 
topology, spaced to effectively span both the massless and degenerate LSP regions. However, 
we found that these benchmarks alone fell far short of our goal. A search optimized only 
for these benchmarks will miss roughly one third of the simplified model parameter space. 
To improve the benchmark list, we found that the most important additions were in the 
simplified model topologies with many b-quarks but without top quarks and leptons, such 
as the Q^^ and Bj^ models. The set of benchmarks found are listed in Appendix [C| Note 
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TABLE I: Search regions that were optimized on the benchmark models. Together they cover 
virtually all models with an efficacy of 1.75 or better, for jO = 1 fb~^, 5 fb~^, and 15 fb~^. 



that the set is not unique, and it may also be possible to create a list with a slightly smaller 
number of benchmarks; however, the list does present a useful solution to covering the 
simplified model parameter space. 

Table |T] presents the search regions optimized on the benchmark models that together 
cover virtually all models with an efficacy of 1.75 or better (i.e. Scrit = 1-75), for C = 1 



fb , 5 fb , or 15 fb . Fig. 10 shows the individual background contributions to the 



optimal search regions found for £ = 5 fb ^. In Fig. 11 , we show for £ = 5 fb ^ the region 



covered by each of the seven signal regions for the Q"^^ and simplified models. Appendix 
[B] shows the expected 95 CL. limits cross section times branching ratio sensitivity for all 
simplified models in Sec. [Tljfrom the search regions in Table |T] that have been optimized on 
the benchmarks in Appendix [c| for C = lfb~^, 5fb~^, and 15fb~^. 

An important observation is that the 36 and SSDL channels should be utilized. Exper- 
imental studies, which currently only go up to 2b and li, could be better optimized by 
considering these channels once data sets with greater than 1 fb~^ begin to be studied. 
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FIG. 10: Individual background contributions to each of the seven search regions in Table |I] for 
C = 5 fb^^. The dominant background is usually tt + nj. The SSDL search region is essentially 
background- free. 



Fig. 11 shows that the 36 and SSDL cuts are by far the most sensitive for top-rich samples. 
The SSDL cut plays an increasingly pivotal role as the size of the data sample increases. The 
SSDL cut is especially useful for reducing background without requiring significant amounts 
of missing energy. The combination of 6-jets and SSDL is crucial in order to achieve optimal 
sensitivity to the simplified model near the degeneracy line where rrig ^ 2mt + m^o. 
Furthermore, even in the case of topologies with mixtures of tops and bottoms, such as the 
Q™ simplified model, SSDL will be one of the dominant discovery channels. The 36 cuts are 
also useful for discovering simplified models with many 6-jets and missing energy, such as 
the Q^i model. 

To further explore the utility of 36 and SSDL channels, we constructed optimized searches 
without either or both, and the results are striking. Without the SSDL channel, the sen- 
sitivity to top-rich samples, especially those with low missing energy, will be significantly 
suboptimal. Searches without 36 channels seem to be even more problematic, as our sensi- 
tivity to virtually all heavy flavor color octet decays is degraded. In fact, we were unable to 
find a set of search regions that are able to cover a majority of the parameter space with an 
efficacy better than 1.75. 

The solutions to the optimization problem depend on the choice of systematic errors 
made. It is natural to ask if a larger systematic error in the 36 region would decrease its 
utility, but we have found that even with e(nbjet)syst = 1, the 36 signal region will still be 
useful for discovering these simplified models. 
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FIG. 11: The regions covered by each of the seven searches in Tablejljfor the (top seven plots) 
and (bottom seven plots) topologies for £ = 5 fb~^. Note that the 3b search is particularly 
effective for these two topologies, and the SSDL search is also important for Q^^ . The black dots 
are benchmark points for these topologies (see Appendix O. 
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Altogether, the inclusion of 3b and SSDL channels leads to a great improvement in efficacy 
for the majority of heavy-flavor simplified model parameter space. Although we cannot 
draw firm conclusions without a full LHC detector study, our results suggest that the LHC's 
sensitivity to new physics will be greatly improved by including these channels. 



VI. DISCUSSION 

This work presented a framework for constructing optimal search strategies sensitive to 
heavy flavor and missing energy signatures at the LHC. We used a set of simplified models, 
with each model or topology parametrized by only a small number of parameters. This 
model space offers a wide range of kinematics and can only be probed with a broad and 
flexible search strategy. A search strategy consists of a list of counting experiments, each 
to be performed in a given search region where a particular set of kinematic and selection 
cuts have been applied to the data. Theorists cannot determine the optimal search strategy 
because we do not have access to realistic LHC detector simulations. However, we can 
sidestep this issue by instead providing a more robust and useful piece of information: a set 
of benchmark models with the property that any search strategy sensitive to all of them will 
also be sensitive to the entire parameter space of heavy flavor simplifled models. 

The benchmarks in Appendix [C] have been designed to span the parameter space of 
simplifled models involving the pair production of color octets or triplets decaying to all 
plausible combinations of tops, 6-quarks, light flavor quarks, and missing energy. Since the 
color octets can decay to any combination of third generation quarks, or a pair of light flavor 
quarks, there are a large number of possible topologies, which differ qualitatively in terms 
of the number and momenta of jets, leptons, and 6-jets, and also the amount of missing 
energy. Thus there is a large parameter space of models, and it is rather non-trivial that a 
small number of signal regions can cover all of these models very effectively. 

The notion of efficacy, S, was used to get a quantitative handle on the optimization. The 
efficacy of a search region applied to a given model is deflned by the ratio of the amount of 
data needed to discover the model to the amount of data needed to discover the model using 
the optimal search region. Our results suggest that it is possible to obtain £ < 1.75 over the 
entire parameter space of heavy flavor simplifled models with a search strategy consisting 
of only 6 or 7 search regions, depending on the integrated luminosity. 

The search regions specifled here are powerful enough to cover many models simultane- 
ously, but one can do better for particular models with a dedicated search. For example, 
to gain sensitivity for t — t + other methods may be required for the 7 TeV LHC to 
separate the signal from the dominant tt background, which is at least six times larger than 



the signal. Fig. [12] shows the and Ht distributions for a 350 GeV stop that decays to a 
50 GeV LSP. Separating the signal from background is quite challenging for this topology 
because the signal peaks at low and is lost amidst the SM backgrounds. Another chal- 
lenging topology is Q^. Fig. 13 shows the and Nj distributions for a 900 GeV gluino 
that decays to a 150 GeV LSP for both the and topologies. Because the highest 
jet-mutliplicity search that we consider is 4+ jets, the cuts we propose for g ^ ti + are 
suboptimal, as more jets would help to reduce background versus signal (see e.g. [311 132]). 

One lesson that emerged from this study is that the character of the optimal search 
strategies will change signiflcantly as larger amounts of data are analyzed. In particular, 
while searches involving 36 jets and same-sign dileptons may not be particularly useful with 
less than 1 fb~^ of data, it seems that with more data these channels will be crucial for 
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EJt (GcV) Ht (GeV) 



FIG. 12: (left panel) and Ht (right panel) distributions for a 300 GeV stop that decays to a 50 
GeV LSP from the topology vs SM backgrounds. The distribution is plotted after requiring 
Nj > 4, TVf, > 1 and Ht > 400 GeV. The Ht distribution is shown after requiring Nj >^,Ni,>l 
and ^T > 200 GeV. 




FIG. 13: f^T (left panel) and Nj (right panel) distributions for a 900 GeV gluino that decays to a 
150 GeV LSP for the G"^^ and topologies. The distribution is plotted after requiring Nj > 4 
and Nb > 3. 



obtaining optimal sensitivity to heavy flavor simplified models involving color octets. This 
conclusion appears to be robust: even if the systematic uncertainties on the 36 backgrounds 
are taken to be as large as 100%, the 36 searches will still be important for fully utilizing 
the LHC data. We look forward to an exciting year as the LHC completes its 7 TeV run, 
taking an order of magnitude more data. 
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of this work, we discovered several similar but distinct studies on heavy flavor and missing 
energy [331435] . 



Appendix A: Expected limits from existing LHC searches: Plots 



In Fig. 



in Sec. IV 



section. 



14, we show the expected 95% CL. limits for £ = 1 fb ^ from the ATLAS studies 



on all simplified models from Sec. Ofor different choices of the production cross 



Appendix B: Expected limits from optimized searches: Plots 



In Figs. 15, 16 , and[T7} we show the expected 95% CL. limits cross section times branching 
ratio sensitivity for the simplified models in Sec. |Tl] from the search regions in Table [T] that 
have been optimized on the benchmarks in Appendix |C[ for C = 1 fb~^, 5fb~^, and 15fb~^. 



Appendix C: Benchmarks 



The list of benchmark models from the space of 12 simplified models is presented in this 
section. The benchmarks from the and Qyy simplified models are presented in Table [IT] 



and Tables III and IV, respectively. The benchmarks from the 7]j and Bj^ topologies are 
shown in Table |Vl 
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FIG. 14: Estimated 95% C.L. contours for £ = Ifb"^ for the cross section times branching ratio 



sensitivity for various simphfied models using the ATLAS jets and searches in Sec. IV The 
soUd, dashed, and dotted Hues give the Hmit for app^xx = 1,3, and 0.3 times c^^xx ' where 
X = g for all plots except for the middle plot in the bottom row (X = t) and for the right plot in 
the bottom row {X = b). 
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FIG. 15: Estimated 95% C.L. contours for the cross section times branching ratio sensitivity for 
various simphfied models using the search regions in Table |l] that have been optimized on the 
benchmarks in Appendix [c| Shown are C, = lfb~^ (left column), 5fb~^ (middle column), and 
15fb~^ (right column). The solid, dashed, and dotted lines correspond to (jpp^gg =1, 3, and 0.3 
times o\lf^;S'^^ ■ Each row is for a different simplified model; from top to bottom, these are G^^ ■, 

qBH qMH qTB 

^ly' V' "^^"^ ■ 



26 




^™5fb-' 









1000 



800 



^ 600 



= 400 



200 



400 600 800 1000 1200 1400 
nig [GeV] 



00 600 800 1000 1200 1400 
m,, [GeV] 



400 600 800 1000 1200 1400 
mj [GeV] 



1000 



800 



^ 600 



400 



200 















— tiA 50 



1000 






800 






600 






400 






200 







400 600 800 1000 1200 1400 
nig [GeV] 



400 600 800 1000 1200 1400 
nu [GeV] 



1000 




/ 200 


800 






600 






400 






200 













400 600 800 1000 1200 1400 
«ig [GeV] 




400 600 800 1000 1200 1400 
mj [GeV] 




00 600 800 1000 1200 1400 
mj [GeV] 



1000 
800 



> 

^ 600 



400 



200 







f^f.ifb-' A 

ummr /, 


















20 



200 400 600 800 1000 1200 1400 
ms [GeV] 




200 400 600 800 1000 1200 1400 
ms [GeV] 



200 400 600 800 1000 1200 1400 
nij [GeV] 



FIG. 16: Estimated 95% C.L. contours for the cross section times branching ratio sensitivity for 

various simphfied models using the search regions in Table |l] that have been optimized on the 

benchmarks in Appendix [c} Shown are C = lfb~^ (left column), 5fb~^ (middle column), and 

15fb~^ (right column). The solid, dashed, and dotted lines correspond to app^gg =1, 3, and 0.3 

times crfl^-?^^. Each row is for a different simplified model; from top to bottom, these are Qj^,, 
PP^yy W 

g^y, gv, and . 
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FIG. 17: Estimated 95% C.L. contours for the cross sectfon times branching ratio sensitivity for 
various simphfied models using the search regions in Table |l] that have been optimized on the 
benchmarks in Appendix [c} Shown are £ = Ifb^^ (left column), 5fb^^ (middle column), and 
15fb~^ (right column). The solid, dashed, and dotted lines correspond to app^xx =li 3, and 0.3 
times o^p^x'x^ 'w^^^re X = g (top two rows), X = i (third row), ot X = b (last row). Each row 
is for a different simplified model; from top to bottom, these are Qf^j, G''^ , T^, and 
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Name 


m~g (GeV) 


m^o (GeV) 


^rft-i (fb) 


(T-^'Li (fb) 


a--l_, (fb) 


<od (fb) 


^B 
^B 
^B 


500 

500 
650 
800 
800 


115 
40 

40 
415 
40 


592 
428 
139 

469 
92 


129 
95 

65 

129 

27 


44 
32 

26 
44 
13 


2310 
2310 

335 
61 
61 


^B 
^B 

g^3 

^B 
^B 
^B 

gB.B 
^B 
gBB 

gB_B 
^B 
gB_B 
^B 
gB.B 


100 
200 
200 
350 
500 
600 
700 
750 
800 
800 
900 


40 
15 
165 
165 
40 
365 
265 
490 
765 
40 
540 


353000 
17800 
3360 
875 

94 
236 

57 
153 
4056 

42 

65 


265000 
11400 
3230 

591 

37 

112 

20 

62 
1840 

11 

23 


226000 
10400 
3210 

373 

24 

70 

11 

41 

1490 
5.2 

13 


21.2x10^ 
625000 
625000 
24200 
2310 

617 

186 

106 

61 

61 

21 


^B 
^B 

gl^ 

^B 

gi^ 

^B 
^B 

^B 
^B 


150 
200 
300 
500 
600 
650 
750 
800 
850 
850 


15 
165 
115 
40 
515 
115 
215 
765 
40 
590 


12900 
39300 
6450 

406 
2590 

129 

90 
3700 

517 

39 


128000 
25700 
4970 

306 
1440 

82 

52 
2750 

351 

19 


115000 
19900 
4300 
278 
939 
67 
41 

2250 
244 
12 


2.86x10^ 
625000 
62100 
2310 

617 

335 

106 

61 

36 

36 



TABLE II: Benchmark models from the pure g^ simplified models. Also shown are the estimated 
cross section reach for £ = 1,5, 15 fb~^ in addition to the NLO production cross section. 
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Name 


m~g (GeV) 


m^o (GeV) 


^ift-i (fb) 


(T-^'Li (fb) 


a--l_, (fb) 


<od (fb) 


^B 
^B 
^B 


500 

500 
650 
800 
800 


115 
40 

40 
415 
40 


239 
175 

88 
152 
66 


146 
100 

29 
59 
17 


92 
63 

14 
37 
8.3 


2310 
2310 

335 
61 
61 


^B 
^B 
^B 


450 
550 
650 
800 
800 


65 
140 
40 
415 
40 


1680 
653 
177 
349 
79 


1320 

470 

102 

234 

39 


1080 
354 
83 
183 
24 


4760 
1170 
335 

61 

61 


B 
^B 
^B 

gv 

^B 

gv 

^B 


200 
200 
500 
800 
800 


165 
40 
40 
765 
40 


25000 
35100 

311 
4120 

58 


17900 
25400 

197 
2960 

29 


13000 
11800 

179 
2510 

17 


625000 
625000 
2310 

61 

61 



TABLE III: Benchmark models from the hybrid g^ simphfied models. Also shown are the estimated 
cross section reach for £ = 1,5, 15 fb~^ in addition to the NLO production cross section. 



Name 


m~g (GeV) 


m^o (GeV) 


^rft-i (fb) 


(7--'Li (fb) 


a--l_, (fb) 


<od (fb) 


gm 

g™ 

gm 


500 
500 
650 
800 
800 


115 

40 
40 
415 

40 


422 

324 
115 
243 
81 


184 

126 

52 

130 

25 


63 

44 
25 
66 
12 


2310 

2310 
335 

61 

61 


gBn 

gm 

gBH 
gBH 
gBH 


300 
400 
600 
800 
800 


45 

220 

170 

595 

45 


1370 
2660 

113 
1160 

55 


1180 
1300 

40 
452 

15 


1010 
619 
25 
240 
6.9 


62100 
10400 

617 

61 

61 


gm 
gm 
gm 
gm 
Qm 


300 
450 
550 
800 
800 


45 

270 

45 

595 

45 


3230 
3190 

150 
1290 

69 


695 
1530 

86 
727 

21 


272 
674 

51 
413 

10 


62100 
4760 
1170 

61 

61 



TABLE IV: Benchmark models from the g^^ simplified models. Also shown are the estimated cross 
section reach for £ = 1, 5, 15 fb^^ in addition to the NLO production cross section. 
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IN dlllc 


I'f^g V j 


If t^O ^VJC V J 


^reach / fU\ 


^reach / fU \ 


-.reach / fU\ 
^15fb-i ^^"^^ 


QCD (r<\ 




250 





15100 


9960 


5980 


180000 




35U 


50 


1970 


1500 


1104 


24200 




500 


200 


536 


349 


289 


2310 




500 


50 


240 


124 


104 


2310 


rs 


650 


350 


321 


178 


144 


335 




650 


50 


96 


49 


32 


335 


^B 
^B 

Bb 

^B 
^B 
^B 
^B 
^B 
^B 


100 
200 
350 

350 
400 
450 
500 
500 
550 
600 




50 
200 

50 
50 
150 
350 
50 

350 


219000 
11200 
2260 

481 
263 
230 
989 
142 
121 
233 


203000 
8620 
1680 

438 
209 
168 
586 

71 

65 
153 


124000 
5370 
1260 

427 
171 
133 
348 

54 

45 
120 


21.2x10^ 
625000 
24200 

24200 

10400 

4760 

2310 

2310 

1170 

617 



TABLE V: Benchmark models from the hybrid and simpUfied models. Also shown are the 
estimated cross section reach for £ = 1,5, 15 fb~^ in addition to the NLO production cross section. 
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